function[agrupacion]=coligamiento(xbin)
[m,n]=size(xbin);
xx=zeros(m,n);
xx(find(xbin==0))=1;
se=[0 1 0;1 1 1;0 1 0];
xx=imclose(xx,se);
[L,num]=bwlabel(xx,8);
valor=[];
for i=1:num
   clear f c
   [f,c]=find(L==i);
   if length(f>30)
       valor=[valor (max(c)-min(c))];
   end
end
agrupacion=mean(valor);